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Abstract. - The dynamics of a single qubit interacting by a sequence of pairwise collisions with 
an environment consisting of just two more qubits is analyzed. Each collision is modeled in 
terms of a random unitary operator with a uniform probability distribution described by the 
uniform Haar measure. We show that the purity of the system qubit as well as the bipartite 
and the tripartite entanglement reach time averaged equilibrium values characterized by large 
instantaneous fluctuations. These equilibrium values are independent of the order of collision 
among the qubits. The relaxation to equilibrium is analyzed also in terms of an ensemble average 
of random collision histories. Such average allows for a quantitative evaluation and interpretation 
of the decay constants. Furthermore a dependence of the transient dynamics on the initial degree of 
entanglement between the environment qubits is shown to exist. Finally the statistical properties 
of bipartite and tripartite entanglement are analyzed. 



k><( " Introduction. — The repeated collision model has 
been recently used in literature to analyze the irreversible 
dynamics of a qubit interacting with a reservoir consisting 
of a large number of environmental qubits. In particular 
processes like thermalization [1] and homogenization [2-5], 
have been analytically investigated. The same model has 
ben used recently also to analyze the dynamics of a qubit 
interacting with a very small environment consisting of 
just two qubits [6]. The interest for such system is due to 
the fact that, at variance with what happens in the case of 
an environment with a large number of degrees of freedom, 
the system dynamics cannot be described by a Markovian 
master equation. Indeed, due to the fact that the sys- 
tem qubit collides repeatedly with the same environment 
qubits, the dynamics is characterized by large fluctuations 
and only when the sequence of collision is random a time 
averaged equilibrium is reached. While in all the above 
mentioned papers the - elastic - collisions have been mod- 
eled by a partial swap unitary operator, in the present 
paper we will analyze the system dynamics in the case in 
which the two-qubit collisions are described by random 



unitary operators [7-12]. 

Random unitary operators have received considerable 
attention in quantum information theory, mainly be- 
cause they find applications in various quantum proto- 
cols [13-16]. Unfortunately the implementation of a ran- 
dom unitary operator acting on the n-qubit Hilbert space 
requires a number of elementary quantum gates that is 
exponential in the number of qubits. On the other hand 
sequences of random two-qubit gates (collisions) generate 
pseudo-random unitary operators which approximate, to 
the desired accuracy, the entanglement properties of true 
n-qubit random states [17-23]. This approach has given 
very good results, showing that pseudo-random states can 
be generated efficiently, that is polynomially in n [17-23]. 

Our choice to describe the pairwise collisions in terms of 
random two qubit unitary operators is motivated by the 
fact that often a precise modelization of the interaction is 
very hard and that, on the other hand, a good descrip- 
tion of the approach to equilibrium can be obtained by 
suitable averages of the quantities of interest, as we will 
describe below. Furthermore such collision model exhibits 
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interesting features ranging from memory effects to the 
efficient entanglement generation between system and en- 
vironment. 

The model. — In order to illustrate the approach of 
our work, we first review the repeated collision model. Let 
us consider a set of iV + 1 qubits, the first of which is the 
system qubit and the remaining N are the reservoir. The 
system-environment interaction is due to pairwise colli- 
sions between the system and a singe reservoir qubit. Af- 
ter t collisions the overall state of the system plus reservoir 
is 

g^E = ■ ■ ■ U.,, U.,, B%UI Ul---Ul (1) 

where qse is the total density operator and the sequence 
ii ■ ■ ■ it specify the order with which the environment 
qubits collide with the system one. As in [6] we have 
concentrated our attention to the case in which the en- 
vironment consists of just two qubits. In the following 
will label the system qubit while 1,2 will label the en- 
vironment qubits. At variance with the previous work, 
however, we have considered here the case in which each 
collision is described by a random unitary operator Ui 
picked up from the uniform Haar measure on the group 
t/(4). In particular in our calculations we have found con- 
venient to parametrize each random unitary matrix Ui , in 
terms of the Hurwitz representation of the unitary group 
f/(4) [12,20,26]. Such different choice of the collision uni- 
tary operators has several consequences. First of all the 
collisions are of course no longer elastic and, regardless 
of the number of environment qubits, homogenization is 
no longer achieved. However, even for a few qubits en- 
vironment, a time averaged equilibrium state with large 
fluctuations is reached regardless of the order with which 
the qubits collide (in the numerical data shown below, wc 
consider random sequences ii,---it). Furthermore such 
equilibrium state is independent of the initial state of the 
qubits. In the following wc will characterize some aspects 
of such approach to a time averaged equilibrium state. 

Random qubit-environment interactions have been re- 
cently considered [24,25], for high dimensional environ- 
ments. However in our work we deal with small environ- 
ments, such that the information acquired by the environ- 
ment on the system can flow back to the system and a 
Markovian description of our model is surely not possible. 
Moreover, there is no weak coupling parameter and the 
state of the environment is significantly affected by the 
interaction at each collision, so that also the Born approx- 
imation does not apply. 

Approach to equilibrium. — We have characterized 
the decoherence of the system qubit in terms of the purity 
V. We remind the reader that the purity is defined as 
V ^ Tr [gl] , where gs — T^^e [qse] is the reduced density 
operator of the system qubit. The purity is a decreasing 
function of the degree of statistical mixture of the qubit 
and takes values in the range < V < 1 where V ^ 1 
corresponds to pure states and V ^ to the completely 



unpolarized mixed state. As already mentioned, due to the 
small number of environment qubits, the instantaneous 
system purity undergoes large fluctuations regardless of 
the initial state of both the system and the environment 
and regardless of the sequence of collisions as shown in 
Fig. 1 (first row). The purity, however, approaches a time 
averaged equilibrium state. To see this wc calculate the 
time averaged purity VrAit) as 

VTA{t)^J^J2^'[Psit')]- (2) 

As shown in Fig. 1 (second row), VxAit) reaches the same 
equilibrium value regardless of the initial entanglement 
of the environment qubits, a natural consequence of the 
random nature of the collisions. 




Fig. 1: From top to bottom: instantaneous, time averaged 
and ensemble averaged purity, for an environment in an initial 
product (left) or maximally entangled (right) state. The time 
evolution of the ensemble averaged purity is in both cases well 
fitted by the exponential curve VEAit) — Vl oc exp(— 0.36f), 
with Vl = I (asymptotic horizontal line). 

Further insight into the approach to equilibrium of the 
system qubit is obtained if rather than the time average 
purity we evaluate the ensemble averaged purity VsAit) 
over the uniform Haar measure. As we shall see, an anal- 
ysis of the time evolution of such averaged quantity allows 
us to determine the time scale for the relaxation to equilib- 
rium. To evaluate VsAit) we have generated a large num- 
ber of sequences of random collisions, each one drawn from 
the uniform measure, and we have averaged the purity over 
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the different histories after t steps. Also such ensemble 
average shows an irreversible behaviour of the qubit dy- 
namics. Regardless of the initial state of the environment 
qubits, VEA{t) decays to the value |, which coincides with 
the average purity Vl predicted by Lubkin [28] for true 
overall (system-environment) random states. Note that 
'pTA{t) and VsAit) tend to the same asymptotic value. 
If fi and v are respectively the dimensions of the Hilbert 
space of the system and of the environment one has 



In our case fi ~ 2 and = 4, and therefore Vl = §■ 
However, as shown in Fig. 1 (third row) the time evolution 
of VsAit) exhibits a clear dependence on the degree of 
entanglement of the initial state of the environment qubits. 
Notably the exponential approach to equilibrium is for the 
first collisions a decreasing (increasing) function of t for an 
environment in an initial product (maximally entangled) 
state. In both cases a numerical fit shows that VEA{t) — 
Vl oc cxp(— At) with A k, 0.36. 

The asymptotic relaxation to equilibrium can be com- 
puted analytically following the approach developed in 
[21-23]. Each pure three-qubit state psE can be 
expanded over products of Pauli matrices: psE = 

Eao,ai,a2 Cao"ia2<° ® ® ^ whcrC Crf dcnotCS a 

Pauli matrix acting on the ith qubit, with ai G {0, x, y, z\ 
and (T° = /. The purity then reads V{i) = X^ao c^oOo(0- 
The purity decay can therefore be obtained from the evolu- 
tion in time of the coefficients c^^qo- ^^"^ '^^'^ random colli- 
sion model the column vector (? of the coefficients „ „ 
evolves according to a Markov chain dynamics of the form 
c2(t 4- 1) = M(?[t) [23]. The matrix M is obtained af- 
ter averaging over the two possible couplings 01 and 02: 
M = +M^^), with M^f acting non trivially (dif- 

ferently from identity) only on the subspace spanned by 
qubits i and j. Averaging over the uniform Haar mea- 
sure on U (4) one can sec that Af^ preserves identity 
((7° ® ® cr^j) and uniformly mixes the other 15 

products ' ® ct"^ [23] . The matrix M has an eigenvalue 
equal to 1 (with multiplicity 2) and all the other eigenval- 
ues smaller than 1. Therefore the asymptotic purity decay 
is determined by the gap A in the Markov chain, namely 
by the second largest eigenvalue 1 — A of the matrix M . 
We have VEA{t) - Vl - (1- A)* ^ exp{[ln(l - A)]t}. 
In our model 1 — A = 0.7 (eigenvalue with multiplicity 2) 
and therefore — ln(l — A) « 0.357, in very good agreement 
with the value A « 0.36 obtained from our fit. 

Entanglement dynamics. — The dynamics of bipar- 
tite and multipartite entanglement between the qubits of 
our model shows interesting features. Such dynamics has 
been conveniently characterized in terms of the concur- 
rence and of the tangles [30,31]. We remind the reader 
that, given the density operator pij of a bipartite system 



of two qubits, the tangle Ti\j is defined as 

'Ti\]ip) = [niax{0,Q;i - Q!2 - as - 04}]^, (4) 

where {ak} [k = 1, ..,4) are the square roots of the eigen- 
values (in non-increasing order) of the non-Hermitian op- 
erator pij = Pij {ay (8) ay) Pij {ay ® ay), ay is the y-Pauli 
operator and p*j is the complex conjugate of pij, in the 
eigenbasis of the az ® a^ operator. The concurrence C is 
defined simply as Cij = yjTi\j. The tangle r^ij, or equiv- 
alently the concurrence Cij can be used to quantify the 
entanglement between the pair of qubits i,j for an arbi- 
trary reduced density operator pij. Furthermore, when 
the overall state of the system is pure, the amount of 
entanglement between qubit i and all the remaining can 
be quantified by the tangle Ti\Y&si ~ 4detpi. The tan- 
gle TQirest between the system qubit and the environment 
conveys the same information as the purity V . Indeed 
it is easy to show that Tg|j.gg^ = 2 — 27^. We have nu- 
merically computed the tangles ro|i, to|2, and ri|2 of the 
two-qubit reduced density matrices and the three-tangle 
Ti\j\k = T%\2k - Ti\j - Tiife, where i,j, k can be any permu- 
tation of 0,1,2 and where the tangle Ti\j^. measures the 
entanglement between the «th qubit and the rest of the 
system, i.e., qubits j, fc. The three-tangle to|i|2 is a mea- 
sure of the purely tripartite entanglement and is invariant 
under permutations of the three qubits [31]. 

The instantaneous dynamics of the tangles Ti\j and of 
the three-tangle To|i|2 is similar to the dynamics shown 
in the Fig. 1 (top) for the purity, that is, these quanti- 
ties are characterized by large instantaneous fluctuations. 
However the time averaged tangles {Ti\j)TA{t), defined in 
analogy with Eq.(2), approach the same limiting value Tp. 

Again we have numerically evaluated the ensemble 
average (with respect to the Haar measure) tangles 
{Ti\j)EA{t)- As shown in Fig. 2 (bottom) the pairwise 
tangles (to|i)£;a(0 and {TQ\2)EA{t) approach exponentially 
the same equilibrium value Tp « 0.367. The numeri- 
cal data in Fig. 2 (bottom) are well fitted by the curves 
{Ti\j)EA{t) — Tp (X exp{—Xijt). When the environment is 
initially separable we obtain Aqi ~ A02 ~ 0.76, A12 « 0.44, 
while for an initially maximally entangled states of the en- 
vironment we have Aqi ~ A02 ~ 0.44, A12 « 0.36. There- 
fore, initially entangled environment qubits limit the rate 
of generation of bipartite entanglement between the qubit 
system and a single environment qubit, even though the 
amount of pairwise entanglement obtained asymptotically 
is always the same. We have furthermore verified also that 
the average tangle {ti\2)ea approaches the same limiting 
value Tp even though the environment qubits do not col- 
lide directly. 

As one would expect from the above discussion, also 
the mutipartite entanglement approaches an equilibrium 
value. This can be seen by the time and ensemble aver- 
ages three-tangle to|i[2 shown in Fig. 3. The approach to 
the equilibrium value tt is exponential also for this quan- 
tity and we can extract the convergence rate from the fit 
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Fig. 2; Time (first row) and ensemble (second and third rows) 
averages for the pairwise tangles ro|i (solid line in the first 
row), ro]2 (dashed line in the first row) and ri[2 (dotted line in 
the first row), for an environment in an initial product (left) or 
maximally entangled (right) state. For the ensemble averaged 
quantities also the fit (Tnj)EA(t) — Tp oc exp(— Aijt) is shown. 
In the second row only the {To\\)EA{t) is shown due to the same 
figure of (ro|2)BA(i) 

(To|i|2)EA(i) - TT (X exp(-Aoi20> with A012 = 0.36 re- 
gardless on whether the initial state of the environment is 
entangled or separable. 

Finally, we note that not only the time and ensemble 
averages of the various tangles converge to the same limit 
values regardless of the initial conditions, but also to a well 
defined limit distribution. In particular the concurrences 
arc distributed in accordance with the distribution shown 
in in Fig. 3 of [27] for random pure 3-qubit states. Fi- 
nally the numerically calculated distributions of the three- 
tangle in our collision model are shown in Fig. 4, for sep- 
arable and for entangled initial environment state. 

Conclusions. — In summary, we have shown that re- 
laxation (in time average) to statistical equilibrium is pos- 
sible for a system of just three qubits undergoing purely 
random unitary evolution, regardless of the order of the 
collisions (in contrast to the case studied in [6] ) . This pro- 
cess is intrinsically irreversible (in contrast to [1,2]) due to 
the random nature of the interactions. The purity of the 
qubit system shows an exponential decay and the decay 
rate has been evaluated both numerically and analytically. 
The limit value is in accordance with [28]. The limit val- 
ues (in a time average sense) of the pairwise qubit tangles 
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Fig. 3: Time (top) and ensemble (bottom) averages for the 
three-tangle ro|i|2 for an environment in an initial product (left) 
or maximally entangled (right) state. For the ensemble aver- 
aged quantities also the fit (To|i|2)BA(t) — tt oc exp(— 0.36t) is 
shown. 




Fig. 4: Statistical distributions of the three-tangle ro|i|2, for 
separable (left) and entangled (right) initial state of the envi- 
ronment. In both cases ro|i|2 — 0.33 ± 0.19. 

are the same for all possible qubit pairs. In contrast to [6], 
residual entanglement (three-tangle) is generated regard- 
less of initial entanglement of the environment qubits. 
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